Quantum statistics of overlapping modes in open resonators 
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We study the quantum dynamics of optical fields in weakly confining resonators with overlapping 
modes. Employing a recently developed quantization scheme involving a discrete set of resonator 
modes and continua of external modes we derive Langevin equations and a master equation for 
the resonator modes. Langevin dynamics and the master equation are proved to be equivalent in 
the Markovian limit. Our open-resonator dynamics may be used as a starting point for a quantum 
theory of random lasers. 
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I. INTRODUCTION 

The interaction of a single resonator mode with an ex- 
ternal optical field leads to damping and noise for the 
intracavity mode. The textbook example is a single os- 
cillator linearly coupled to a continuum of harmonic os- 
cillators 0,0, HIl- Less understood and in fact object of 
ongoing debate |E IE 0, IE 0, E3, El E3, are the damp- 
ing and noise properties of multimode fields in resonators. 
Multimode fields show an excess noise or Petermann fac- 
tor |E 0, 0| • In laser systems this factor gives rise to 
a peculiar enhancement of the laser linewidth above the 
fundamental Schawlow-Townes value; this enhancement 
was measured in recent experiments on unstable laser 
cavities 0, . Excess noise and the enhancement of 
the laser linewidth may be attributed to the nonorthogo- 
nality and the spectral overlap of the cavity eigenmodes 
in the presence of the coupling to the external field. 

While the phenomenon of excess noise has been known 
for more than 20 years, there is still no complete descrip- 
tion of the quantum statistics of overlapping modes. A 
few example systems were discussed H,l3| from a quan- 
tum mechanical point of view. Recently several authors 
|lOl ITlL Il2l Il3| | proposed quantum Langevin and mas- 
ter equations for multimode fields; however, the status 
of these equations remained unclear as they were not 
derived from rigorously quantized electromagnetic fields. 
This is in contrast to the quantum properties of a single- 
mode field (linearly coupled to an external heat bath) 
which are known |iJ. y. |jj| for arbitrary damping strength 
and arbitrary heat-bath temperature. 

The goal of the present paper is to derive and clarify 
the status of stochastic equations for the field dynamics 
in resonators with overlapping modes. The experimen- 
tal motivation for our work derives both from realiza- 
tions |u| ^| of unstable laser cavities and from recent 
experiments of highly disordered dielectrics which form 
mirrorless so-called random lasers 0, ■ We address 
the field dynamics both within the Heisenberg picture 
(in terms of quantum Langevin equations) and within 
the Schrodinger picture (employing a master equation for 
the reduced density matrix of the cavity modes). We go 
beyond previous work in the following respects: (i) We 
derive the field dynamics starting from rigorously quan- 



tized electromagnetic fields. In particular, no restriction 
of the dimensionality and the vector character of the field 
strengths are indulged in. Keeping only resonant terms in 
the field Hamiltonian and adopting a Markov approxima- 
tion for the memory of the "bath" (i.e., the electromag- 
netic field outside the resonator) , we obtain the equations 
of motion for the intracavity field, including expressions 
for all damping and noise forces, (ii) Our derivation clar- 
ifies the status of the resulting stochastic equations. In 
particular, we show that they correctly describe sepa- 
rated as well as spectrally overlapping resonances. Cor- 
rections would only be important if the resonance widths 
were to become comparable to the resonance frequencies 
themselves, a regime not encountered in optical or even 
microwave resonators, (iii) We derive a representation of 
our master equation in terms of non-orthogonal modes 
and compare our result with master equations proposed 
earlier [lCl fllL Il3| . This allows us to specify the physical 
conditions under which those previously proposed master 
equations hold. 



II. FIELD QUANTIZATION FOR 
OPEN-RESONATOR GEOMETRIES 

There is nothing to add to the familiar canonical quan- 
tization of the electromagnetic field 0,0- The field com- 
ponents may be expanded in terms of any complete set 
of basis functions. As a matter of convenience one usu- 
ally employs eigenmodes pertinent to the given geometry 
and respecting the physical boundary conditions. Such 
expansions are at issue here. 

In the presence of a more or less open resonator, one of- 
ten wants to distinguish between "inside" and "outside" 
even though any opening, i.e., a hole in the material walls 
of the resonator gives to the latter concepts an element of 
arbitrariness. Any choice of a fictitious surface covering 
the hole yields its own inside/outside separation. More- 
over, different boundary conditions may be imposed at 
the chosen separating surface. Nevertheless, each such 
surface and boundary condition entail eigenmodes allow- 
ing us to represent the electromagnetic field almost ev- 
erywhere, with the qualifier "almost" reminding us of the 
fact that the expansion cannot be expected to converge 



2 



pointwise everywhere, and in particular not on an arbi- 
trarily chosen boundary. 

The freedom in choosing the separating surface and 
the boundary condition thereon may be used 
to define a discrete set of inside modes which vanish out- 
side; and similarly a continuum (or even several continua 
distinguished by a "channel" index) of outside modes 
which vanish inside and fulfill scattering-type boundary 
conditions at infinity; a scattering condition specifies a 
single channel for a wave coming in from infinity and 
entails different amplitudes for the partial waves going 
out through the various channels. Canonical quantiza- 
tion then amounts to representing the coefficients of the 
mode expansion of the electromagnetic field by creation 
and annihilation operators {a A ,a^} of photons of the A- 
th inside mode and likewise {b m (u>), b^co)} of photons of 
the (m, w)th outside mode, where m labels channels (in- 
cluding polarization) and u> distinguishes modes within 
the channel continuum. The bosonic commutators read 

[a\,a\,] = 5 X ,\> , (la) 
[b m (oj), b^iui')] = S mm 'S(oj - uj') ; (lb) 

moreover, inside operators commute with outside ones. 

It was shown in two previous papers [20j, |21| that the 
Hamiltonian of the electromagnetic field can be rigor- 
ously expressed as the following bilinear form in the fore- 
going creation and annihilation operators, 

A m J 

J duj{W Xm {Lj)a f x b m (uj) 

A,m 

+Vx m (u)axh m (u) + R.c]. (2) 

Here, the inside mode frequencies lj\ as well as the inside- 
outside coupling amplitudes W\ m , V\m reflect the choices 
made for the separating surface. The latter amplitudes 
are integrals, taken over the separating surface, of prod- 
ucts of inside (or outside) mode functions with deriva- 
tives of outside (inside) mode functions. Specifically, if 
the von-Neumann boundary condition 

n x [V x u A ] \ gl = (3) 

was imposed for the inside modes u\ and the Dirichlet 
boundary condition 

nxv m (u)| s/ = (4) 

for the outside modes v m (oj) at the boundary dl (with 
normal vector n), the coupling amplitudes are given by 

W Am H = -^= f d 2 ru*(r) • {n x [V x v TO (w, r)]},(5) 
2V w aw J dI 

V Ato (w) = — C -== /d 2 ru A (r)-{nx [V X v m (u, r)]}.(6) 
2 v / w A w J dI 



If the physical boundary conditions on material surfaces 
entail time reversal invariance, the amplitudes W, V can 
be chosen real and become identical, W = V. 

The derivation of the Hamiltonian takes full account of 
three spatial dimensions as well as of the vector character 
of the electric and magnetic fields E(x, t) and B(x, t); in 
brief, the Heisenberg equations of motion generated by 
that H arc equivalent to the quantized Maxwell equa- 
tions. Moreover, by composing the inside (outside) mode 
functions with coefficients {a A , a\} ({b m (oj), b^ui)}) one 
gets the fields E(x, t) and B(x,t) inside (outside); trun- 
cated such compositions will be robust against small 
changes of the separating surface unless an unreasonable 
choice was made to begin with. 

Hamiltonians of the above form have been used 0, in 
quantum optics for high-quality resonators; there, inside 
excitations are but slowly dissipated to the outside, with 
life times long not only compared to typical internal pe- 
riods 2ir/u>\ but even to beat periods 2tt/Auj, where Aw 
denotes the modulus of the typical frequency spacing of 
neighboring internal modes. In such applications the an- 
tiresonant terms V\ m a\b m + Vt atfrt, can be neglected. 
However, the derivation of Refs. |20ll21| also secures va- 
lidity of Hamiltonian J5| for very much open resonators 
where the outside field causes the inside resonances to 
overlap, the main case of interest in the present paper. 
The Hamiltonian @ even remains valid in the extreme 
case of overdamped inside excitations, where mode fre- 
quencies uj\ are overwhelmed by large escape rates ft; it 
is only in that extreme situation, which appears as not 
of interest in optics, that the antiresonant terms would 
be important. 



III. LANGEVIN EQUATIONS FOR 
OVERLAPPING RESONANCES 

While Hamiltonian @ describes the coupled dynam- 
ics of the inside and outside fields, its principle appli- 
cation is to separate descriptions of those two subsys- 
tems. As already indicated above, the continua of outside 
modes tend to act as a "bath" damping the discrete in- 
side modes. The effectively irreversible dynamics of the 
inside modes becomes manifest when the bath degrees 
of freedom are eliminated. The Heisenberg equations of 
motion of the inside amplitudes a A , a\ then take the 
form of Langevin equations, in which the outside ampli- 
tudes b m (ui), b^oj) enter only with their initial values 
in inhomogeneities, the Langevin noise forces. Due to 
the bilinear form of Hamiltonian the Langevin equa- 
tions can be determined rigorously by diagonalization of 
H , i.e., without resorting to perturbation expansions (for 
a single "system" oscillator the diagonalization was per- 
formed in, e.g., Refs. [z| and 3). 

The Langevin equations to be noted in the present sec- 
tion are simplified in three respects. First, we bar all an- 
tiresonant terms, assuming the absence of overdamping. 
Second, we confine the discussion to a Markovian situa- 



3 



tion valid for times larger than "bath correlation times" 
T bath (like the thermal time r^ th = fi/fc^T); the exis- 
tence of a time scale separation is thus assumed, such that 
all inside lifetimes are much in excess of Tbath- Third, to 
save space we do not bother to pedantically write out the 
so-called frequency-shift terms which are rarely needed 
in practice. The limit in question yields the Langevin 
equations |21| 

oaCO = -n>\a\(t) - TA/^mW + . ( 7 ) 

wherein the damping matrix 7> M and the noise force Fx 
are given in terms of the inside-outside coupling ampli- 
tudes W\m as 

7A M = t(WH% , (8a) 

Fx{t) = -ifdLJe-^-^Y^^rnbrn^M-i^) 

m 

Note that the damping matrix 7 is non-negative and 
Hermitian; the neglected frequency-shift terms would 
amount to an anti-Hermitian addition to 7. Under con- 
ditions of time reversal invariance the matrix 7 is even 
real symmetric. The Markovian limit mentioned allows 
us to drop the frequency dependence of the coupling am- 
plitudes W over the spectral range in consideration and 
to restrict the time span passed since the initial moment 
to as t — to 3> Tbath- In that limit the force has a white 
spectrum according to (F\(t)) = 0, and the noise 

{Fl(t)F^)) = 2^xrhhS(t-t'), (9a) 
(F x (t)Fl(t>)) = 2 7v (l + n th ) ,$(*-*')• (9b) 

The second two-time correlation functions follows from 
the first and the commutation relations (|lb|) . The re- 
maining second moments vanish, (F\(t)F^(t')) = = 
{F^(t)F^(t')), and higher-order moments follow from the 
ones of orders 1 and 2 according to Gaussian statistics. 
The thermal number of photons n t h = [cxp(?icj/fcT) — 
appearing in the second moments must be taken as 
frequency independent throughout the spectral range of 
inside frequencies oj\ under consideration. In the limit 
kT <C hu) we recover the Langevin equations of Bardroff 
and Stenholm [Io| . 

The appearance of a nondiagonal damping matrix 
7 = 7r>VyV^ signals that the Langevin equation may 
legitimately be applied to the case of overlapping res- 
onances, in which typical matrix elements 7a m are larger 
than a typical nearest-neighbor spacing Aw of frequen- 
cies. However, all elements of 7 must be smaller in magni- 
tude than the frequencies u\ themselves as antiresonant 
terms were dropped in the derivation of Eq. (J7J. Con- 
versely, we could specialize to the weak-coupling regime 
I7A/J *C Aw. Then, lowest-order perturbation theory 
simply amounts to dropping the off-diagonal elements of 
the damping matrix, 7^ — > for A 7^ fx, whereupon the 
Langevin equation J7J simplifies so as to describe a set 
of mutually independent damped harmonic oscillators. 



For some applications it will be helpful to rewrite the 
above Langevin equation with the non-Hcrmitian matrix 



Hx^ = ftujxSxf, - ihjx^ (10) 

diagonalized. The eigenvalues of 7i will then represent 
"true" resonances of the cavity. As a "penalty" for that 
change of representation one would have to work with 
non-standard commutation relations for the operators 
connected with the eigenvectors of Ti. (see Sec. IVIf) . 

A final remark on the dynamics of the outside field is 
in order here. The time scale separation assumed in de- 
riving the above Langevin equations is effective outside 
as well as inside as the outside field evolves in adiabatic 
equilibrium with the inside field for times t — to 2> Tbath- 
For the quantitative treatment by the quantum optical 
input-output formalism or, equivalently, scattering the- 
ory we refer the interested reader to 0, d |2l| ■ 



IV. MASTER EQUATION FOR OVERLAPPING 
RESONANCES 

We now come to the central section of the present pa- 
per. Switching to the Schrodinger picture we present the 
master equation for the reduced density operator p(t) of 
the inside field equivalent to the Langevin equation of 
the preceding section, 

P =-i^2ux[a[ax, p] + (1 + n th ) ^ 7a m |K,pa|] 

A Xfi 

+[<W>4]} +rhh^2~/Xv |[4,pa M ] + [aj^ajjj . (11) 

A/i 

This equation generalizes the familiar quantum optical 
master equation for a single damped harmonic oscilla- 
tor to many oscillators coupled by the (off-diagonal ele- 
ments of the) damping matrix 7^. The latter coupling is 
important when the damping is strong enough to cause 
spectral overlap of modes. The first double sum, pro- 
portional to 1 + nth, describes spontaneous and induced 
emission of photons towards the outside while the sec- 
ond double sum, proportional only to n t h, describes ab- 
sorption from the outside; that interpretation is easily 
checked by employing the Fock representation, i.e. the 
representation in terms of eigenstates of the photon num- 
ber operators a\ax- 

Systematic and stochastic forces are not as clearly sep- 
arated here as in the Langevin equation. In order to bring 
about such distinction here as well, and to prepare for the 
proof of equivalence of Eqs. Q and ltTT|) we may imagine 
the density operator p(t) at time t antinormally ordered 
in the annihilation and creation operators (all a's to the 
left of all cr's); further, we rearrange the commutators 
in the right-hand side of the master equation such 
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that p becomes anti-normally ordered provided p(t) is, 



P = 



-\^2,u>\{[a\, pa\] - [a x p, a{}} 

A 



+2nth^2l\ t ,[[a^,p},a{] 

Xfi. 



(12) 



The latter form of the master equation preserves anti- 
normal ordering of p(t) at all times. Moreover, we have 
now separated reversible drift terms (proportional to the 
frequencies lu\), irreversible drift terms (oc 7a^) not in- 
volving the thermal number of quanta n t h, and noise gen- 
erated diffusion terms (oc n t h); the latter interpretation 
will become obvious in the following section where a rep- 
resentation based on coherent states will be employed. 



V. EQUIVALENCE OF LANGEVIN AND 
MASTER EQUATIONS 

There are various ways of demonstrating the equiva- 
lence of the Langevin equations J7J and the master equa- 
tion and (|12J) . A somewhat indirect (and, in fact, the 
most laborious) way would be to start from the Liouvillc- 
von Neumann equation p in eout(0 = -i[H, Pm®out(t)]/h 
for the density operator of the full electromagnetic field 
and stick to the Schrodinger picture in eliminating the 
outside field, employing the same approximations as 
in deriving the Langevin equations. A standard, and 
more economical, procedure is to show that Eq. J7J and 
(|12[1 entail the same evolution equations for all mean 
values (IIa^a)™* ri AI ( a M)"")(*) with integer exponents 
{to a, n^}. We here prefer to employ a less familiar but 
particularly elegant method. 

If we consistently stick to antinormal ordering of pit") 
we may write the commutators in the master equation 
(|12|l as differential operators as [a, (•)] — * (d/da^)(-) and 
[(•)> a 4 — » (d/da)(-). We may then just as well de- 
grade all creation and annihilation operators to complex 
c-number variables and the density operator to a real 
function of those variables, 



«A 



a* x , p(t)^P({a,a*},t). (13) 



The function P in question is the familiar Glauber- 
Sudarshan quasiprobability [lj which has as its moments 
expectation values of normally ordered observables, 



n( a ir A IlK) n ")w= (i4) 

A A 1 

J J [] da x da* x (a* x ) m > (a A ) n * 1 P({a, a*}, t). 



The master equation l|12f> then becomes the Fokker- 
Planck equation (with the shorthand d/da\ = d x , 



8/dal = d* x ) 

P= [ -iJ2"x(d* x a* x -d x a x ) 

A 

+ ^2 7A M (^a* + df,a x + 2n th d^dl 

X/.i 



(15) 



P. 



At this point the interpretation of the various terms in 
the master equation i|12|) given earlier becomes obvious. 

In the same vein we may degrade the Langevin equa- 
tion to the c-number equation 



OA 



(t) = -~^H x ^(t) + <p x (t), 



(16) 



where tp x (t) is the c-number representative of the oper- 
ator valued random force F x (t). Gaussian statistics and 
zero mean are imparted to <p(t). Moreover, since observ- 
ables are consistently to be taken as normally ordered we 
must identify second moments as 



(<p* x (t)Mt')) 



(Fl(t)F,(t')) 



2n^ lllX 5{t-t')- (17) 



while holding on to ((fi X {t)) — 0, {ip^^ip^t')) = and 
Gaussian factorization for the higher-order moments. 

At this point the equivalence proof can be carried out 
by deducing the Fokker-Planck equation Ijl5(l from the 
c-number Langevin equation (|16fl . and to that task we 
now turn. In order to make the argument as transpar- 
ent as possible we drop mode indices; the initial time at 
which statistical independence of the inside and outside 
fields is assumed is taken as to — 0. We start by con- 
sidering a fixed realization of the time-dependent "ex- 
ternal forces" (p(t) , ip* (t) . The solution of the Langevin 
equation is then a functional of the temporal course 
taken by the external force between the initial and the 
current time. Correspondingly, we introduce a density 
W(a, a*, t\{p(t), <p*(t)}) functionally depending on the 
realization of the external field under consideration; at 
the initial time t = the new density coincides with the 
Glauber-Sudarshan P-function, W(0) = P(0). The den- 
sity W obeys the evolution equation 



W{t) = (L + l(t))W(t), 

Lq = —iui(d*a* — da)+-f(d*a* 

i(t) = -(0V(t) + ap(t)). 



(18) 



da) 



Clearly, we encounter a first-order partial differential 
equation since the random character of the external force 
is not yet invoked. The time independent part Lq of the 
generator accounts for oscillation and damping, while l(t) 
includes the time dependent forces f(t), <p*(t). The for- 
mal solution reads W(t) — (exp J Q ds(L n + i(s))) ,W(0), 
where (•)+ denotes positive time ordering. Temporarily 
employing an "interaction picture" we split off the mo- 
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tion generated by Lq as 
W[i) = e Lot W(t), 
W(t) = [exp / dsl{s) P(0), 



(19) 



Z(t) = e- Lot l{t)e Lot 



[p*{t)d*{t) + <p{t)d{t)] 



Now we may take the average over the Gaussian ensem- 
ble of realizations of (pit) , <p* (t) . The average of the in- 
dependent density W(t) is just the Glauber-Sudarshan 
P-function, (W(a, a*, t\{<p(t), ip* {£)})) = P(a,a*,t). In- 
voking the above first and second moments we get 



P(t) 



exp 



dsl(s) 



P(0) 



exp 



exp 



ds / ds'<<n*(s)^(s'))9*(s)9(s') 



2n th7 / d S d*(s)d(s) 
o 



P(0). 



P(0) 
(20) 



Upon differentiating with respect to time we have Pit) = 
2n th jd*(t)d(t)P(t), and going back to P(t) = e Lot P(t), 
we arrive at the Fokker-Planck equation l|15f) and have 
thus arrived at our goal. 

We should add that a mathematically more satisfac- 
tory variant of the foregoing considerations would result 
from employing the so-called Ito calculus [| . 



VI. STATIONARY SOLUTION OF THE 
MASTER EQUATION 

The general time-dependent solution of the Fokker- 
Planck equation l|12|) is, due to the linearity of the drift 
coefficients and the constancy of the diffusion matrix, 
easy to construct j2^| . We are in fact facing a stochastic 
process of the Ornstein-Uhlenbeck type. The Gaussian 
distribution of the noise together with the linear evo- 
lution equation (|16f) imply that the stochastic variables 
must be Gaussian distributed. In particular, the station- 
ary P function is immediately checked to be 



P({a, a*}) = TT — !— exp(-a* x ax/nth)- 

J- J- 7TT1 



(21) 



The dissipative coupling of the system modes is no longer 
visible in the stationary state; rather, we encounter the 
thermal equilibrium state one would also find in the ab- 
sence of spectral overlap. 



VII. STATUS OF LANGEVIN AND MASTER 
EQUATIONS AND RELATED LITERATURE 

Multimode fields are more complex than single-mode 
fields due the additional frequency scale set by the mean 



frequency spacing Alo of internal modes. For resonators 
with overlapping modes the mode spectral broadening is 
comparable to Alu. Langevin or master equations may 
then only be used if they provide a nonperturbative de- 
scription of damping and noise (in the sense that the 
mode decay rates k may exceed the mode spacing Aw). 

To show that the Langevin Q) and the master equa- 
tions l|ll[l and l|12f> are non-perturbative in the above- 
mentioned sense, we now drop the rotating-wave and 
Markov approximation and address the exact field dy- 
namics. From the field Hamiltonian J5J one obtains the 
following set of Langevin equations 



a\(t) 



-iu\ax(t) - V / dt' [Tx^t - t')a^t') 
\-X Xl .(t-t')al(t'))+fx(t). (22) 



Non-Markovian effects and corrections to isolated-mode 
behavior are encoded in the autocorrelation functions 



iv(i-o = Yl 



dw 



-^-^Wx m (u)w; m (u,) 



+-(*- t V; m (c)^ m (c)l , (23) 



du> 



^w Xm (ui)v: m {u) 



-e +Mt - 4 V; m MiF; m M] . (24) 

The noise term takes the form 



fx(t) 



du> 



- iu{t - to) Wx m (u)b m (u,to) 



+Mt-t ) v * m{uj)b l M 



(25) 



Due to the bilinear form of the field Hamiltonian the 
exact Langevin equations are linear in the system and 
bath operators. In principle, their solution can be found 
by Laplace transformation. The details of the exact so- 
lution will depend both on the spectral strength of the 
bath and the frequency dependent coupling amplitudes 
Wx m (u) and V\ m {w). 

The result (I22|l - H25|) holds for arbitrary damping. 
Here, our main focus is on quantum optical systems 
where the oscillation frequencies Co much exceed the typ- 
ical damping rates k. Then our result can be used to 
compute corrections to the Markovian dynamics of Sees. 
IIIIhTVl We substitute the ansatz a\(t) ~ exp(— i[a> — in]t) 
with the damping rate k into Eq. 12211 ; here ui includes a 
frequency shift 8lo. One obtains four terms proportional 
to <ip or , out of which three terms are strongly oscillat- 
ing: the contribution cx ( W*) has the oscillatory 
integrand <~ 



2iQt' r e +2iut' 



V*W oscillates 



), respectively, and the term 
with respect to ax(t). To 
leading order in the ratio |Jf2|/tD = \SQ — ik\/u> only the 
term cx WW* must be kept. The corrections of higher 
order in |5fi|/w are very small for the systems of interest 
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in quantum optics. We note that there are no correc- 
tions of order |#0|/Aw; thus our field dynamics correctly 
describes the regime of overlapping modes. 

So far we described the field dynamics in terms of 
eigenmodes of a conveniently chosen closed system. The 
associated mode operators {a a, a\} obey canonical com- 
mutation relations. The master equation may equiv- 
alently be expressed in terms of a biorthogonal set of 
modes representing the resonances of the open cavity. It 
turns out that the operators associated with the reso- 
nances obey nonstandard commutation rules. We focus 
on the case kT <C hui and write Eq. in the form 

P = -ji H ^p- pHl s ) + 2^ 7 V a„p4, (26) 
Y, H ^ a W- (27) 

A/; 

The non-Hermitian matrix TL was introduced in Eq. (|10|l . 
The complex eigenvalues of TL represent the resonances of 
cavity in the presence of the coupling to the bath. Gen- 
erally, the eigenvalues are nondegenerate, and TL can be 
diagonalized by a similarity transformation TL — TQT -1 . 
The diagonal matrix Q comprises the eigenvalues of TL 
while T is a general matrix with complex entries. In 
terms of two novel sets of operators 

dn = $>«A°A, (28) 

A 

4 = 5>tT Anj (29) 

A 



the effective Hamiltonian takes the diagonal form 

H cS = ttneUn ■ (30) 

a 

This result resembles the diagonal field Hamiltonian of 
closed resonators. However, the peculiar properties of 
the open resonator dynamics are encoded in nonstandard 
commutation relations of the new operators: From Eqs. 
H28JI and l)29|l and the canonical commutation rules for 
the operators {o,A,a A } one finds, e.g., [c? m ,ejj = 5 m ,n 
and [e m ,ejj = A mn with A = T^T. Substitution of Eq. 
Q28JI into Eq. (|26(1 yields the master equation 

P = -^(H cB p - pHl s ) + ^ ^2 A nm (fl m ~ Q^dmpdl . 

nm 

This is the master equation proposed by Lamprecht and 
Ritsch 13j. Our derivation provides a microscopic basis 
for that equation [24j and proves the equivalence to the 
master equation (|llfl in the limit kT <C huj. 
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